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'?IND TUNNEL SIMULATION OF MARTIAN SAND STORMS 


Objective 

To conduct laboratory simulations of martian windblown processes using 
the Martian Surface Wind Tunnel at NASA Ames Research Center. 

Background 

The Martian Surface Wind Tunnel at NASA Ames Research Center, designed 
and operated by the Principal Investigator, was utilized to carry out a 
broad research program to understand the physics and geological relation- 
ships of particles driven by the wind under near-martian conditions. 

This effort was part of a general program of research on aeolian activity 
as a planetary process funded by the Office of Planetary Geology through 
the Space Science Division at NASA Ames. 

Many physical and chemical process ss modify planetary surfaces, On 
planets having an atmosphere, wind plays an important role in shaping the 
surface through the redistribution of fine grained material; bedrock may 
be eroded, and the particles transported for deposition in other areas. 

This process is readily observable on Earth today, and is evident in the 
geological record. 

Even in its tenuous atmosphere it is now clearly demonstrated that 
aeolian processes are important on Mars. Mariner 9 results show abundant 
wind- related landforms including dune fields and yardangs. Viking lander 
images show drifts of windblown particles; Viking Orbiter pictures add 
to the inventory of aeolian features. On Venus, measurements made by the 
Soviet Vevera spacecraft show winds that may be capable of transporting 


materials and suggest that aeolian processes may operate on that planet as 
well. In addition, one of the satellites of Saturn, Titan, appears to 
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have an atmosphere and atmospheric surface pressures comparable to Mars, 
which would again suggest the possibility of aeolian activity. 

In order to understand the present surface dynamics and the geologi- 
cal evolution of surfaces for planets subjected to wind, it is important 
to know the nature of aeolian processes. Factors such as threshold speeds 
(minimum wind needed to move particles), rates of erosion, trajectories 
of windblown particles, and flow fields over various landforms are impor- 
tant aspects of the problem. 

Threshold Experiments 

The Martian Surface Wind Tunnel was successfully operated during the 
grant period to carry out a variety of studies. Rod Leach, assisted by 
two students from Santa Clara (Mark Plummer and Mary Niemiller) operated 
the facility. Bruce White and his students from the University of Califor- 
nia at Davis were involved with the particle threshold experiments using 
a carbon dioxide atmosphere. Experimental results were combined with 
theory by J. Iversen and J. Pollack to account for temperature difference, 
and a series of curves derived to show minimum wind speeds needed for 
particle motion to occur on Mars. Results were published in early 1980 
(see attached). 

Erosion Rate Experiments 

Preliminary experiments were carried out to determine rates of wind 
erosion of geological materials. Samples of rocks and minerals were bom- 
barded with grains under simulated windblown conditions under a range of 
atmospheric pressures and wind speeds following the procedure outlined 


below: 
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1. Rock samples were cut 5 mm thick, sized to the MARSWIT 
Erosion Facility. 

2. Slabs were numbered for identification and weighed with 
a precision balance. 

3. Samples were bombarded with various particles, crushed 
and sized to an average diameter of 160 micrometers. 

T\vo types of runs were conducted. Case A runs were at 
1 atm. pressure with "Earth air" composition and at 
"Earth" threshold velocities, Case B runs were at 6 mb 
pressure, "martian air" composition, and at martian 
threshold velocities. 

4. Samples were then weighed to determine the loss of weight 
due to erosion, and an erosion rate derived for the two 
cases. 

5. Results were reported at the Lunar and Planetary Science 
Conference (see attached). 


Electrostatic Results 

Effects of electrostatics on particles in the aeolian environment 
were studied and a model derived for windblown materials on Mars, published 
in the Journal of Geophysical Research (see attached). 
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Abstract, Viking observations suggest abundant silt and clay 
particles on Mars, It is proposed that some of these particles 
agglomerate to form sand size aggregates that are redeposited as 
sandlike features such as drifts and dunes. Although the binding 
for the aggregates could include salt cementation or other 
mechanisms, electrostatic bonding is considered to be a primary 
force holding the aggregates together* Various laboratory exper- 
iments conducted since the 1 9th century, and as reported here 
for simulated Martian conditions, show that both the magnitude 
and sign of electrical charges on windblown particles are 
functions of particle velocity, shape and composition, atmu* 
spheric pressure, atmospheric composition, and other factors. 
Electrical charges have been measured forsaltating particles in 
the wind tunnel and in the field, on the surfaces of saml dunes, 
and within dust clouds on earth, Similar, and perhaps even 
greater, charges are proposed to occur on Mars, which could 
form aggregates of silt and day size particles, Electrification is 
proposed to occur within Martian dust douds, generating silt- 
clay aggregates which would settle to thf; surface where they 
may be deposited in the form of sandlike structures, By analog, 
silt-clay dunes are known in many parts of the earth where silt- 
clay aggregates were transported by saltation and deposited as 
‘sand, 1 In these structures the binding forces were later 
destroyed, and the particles reassumed the physical properties 
of silt and clay, but the sandlike bedding structure within the 
Mimes* was preserved. The bedding observed in drifts at the 
Viking landing site is suggested to result from a similar process 
involving silt-clay aggregates on Mars. 

Introduction 

Aeolian (wind related) processes are recognized as important 
in modifying the surface of Mars, as evidenced by observations 
of active dust storms and various erosional and depositional 
features commonly attributed to wind action. However, after 
nearly I decade of intensive Mars studies, there are still many 
fundamental questions left unanswered about the aeolian 
regime of Mars, This report considers the existence of aeolian 
silt-clay aggregates on Mars and presents a model for their 
formation. 

One of the more intriguing results From the Viking landers is 
the apparent lack of sand size particles at the two landing sites. 
The fine-grained material seen on images (Figure 1) is estimated 
to be considerably less than 100 m m in diameter [H. J, Moore et 
al,, 1977] , consistent with estimates of sizes of dust particles in 
the atmosphere [Pollack et ah, 1977], Yet, observations of 
bedding structures in the drifts (Figure 1 ) and of numerous and 
extensive dune fields on Mars argue for the existence of sand 
size particles [Cutts and Smith, 1973], because dunes are 
known to form only from sand size material. Dunes are built 
from the deposition of sal taring particles (and particles in 
tractioi ), and on earth the particle most commonly transported 
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by tation is sand size IBagnoId, 1941], Although ‘sand 1 
is de» led as grains about 60-2000 pm in diameter, most dune 
sand averages 150-300 pm in diameter. This is partly a reflec- 
tion of the wind speeds needed to initiate movement of grains 
(threshold wind speeds); as Bagnold [1941] demonstrated, on 
earth the size most easily moved by the lowest wind speed Is 
aboat 80 pm. Experiments performed under simulated Martian 
conditions [Greeley et ah, 1976] show that the corresponding 
grain size for movement by winds on Mars is abort 160 pm in 
diameter. However, this size particle requires an order of magni- 
tude stronger winds, primarily because of the less dense atmo- 
sphere on Mars, Therefore one may still expect sand size 
particles to be the most active grain size moved by wind and to 
be moved by saltation, but the saltation trajectory would be 
longer, since the winds would be stronger (as demonstrated by 
White et ah [1976] ). TIius dunes on Mars also should be com- 
posed of sand size particles because the basic mechanics of dime 
formation are the same* even taking into account differences in 
the Martian environment. 

The question of sand on Mars is further complicated by 
consideration of potential sources. Smalley and Krinsley 
[1979] point out that most sand on earth is derived from 
granite, a source that is probably lacking on Mars, It is con- 
ceivable that impact cratering of basalts and other rocks on 
Mars would produce u wide range of fragment sizes, including 
sand. Sagan et ah [1977], however, have suggested that under 
the high winds needed to move particles in the Martian atmo- 
sphere, the grains would have a short lifetime, leading them to 
propose the term ‘kamikaze’ grains. In effect, once the sand is 
set into motion by the high winds on Mars, the grains would 
smash into bedrock, other grains, and other objects and be 
fragmented into tiny (e,g. t silt and clay-sized) particles. 

On the basis of laboratory experiments, possible terrestrial 
analogs, and Viking results, a model is proposed in which silt 
and clay-sized particles agglomerate to form sand size grains on 
Mars. 

La b ora to ry li xpe rime n ts 

A series of laboratory experiments was recently undertaken 
to study rates of aeolian erosion of various rock samples under 
simulated Martian conditions using both wind tunnel facilities 
and specially designed apparatus [Greeley and Williams, 1979]. 

These experiments were designed not only to examine 
erosion of individual target particles and rock surfaces but also 
to study effects on the particles doing the abrasion. As Sagan 
et al. [1977] predicted, sand size particles have very short 
lifetimes, given the high velocities and nearly ‘airless’ aeolian 
environment of Mars* Figure 2 shows experiment results. At a 
particle velocity slightly above Martian threshold wind speed, 
sand size material is reduced to particles less than 20 jum in 
diameter (50% by weight) in about 20 min. Figure 2 also shows, 
however, that with continued activity the 50% particle size 
increases after about 60 min* Examination of these particles 
with a scanning electron microscope (5EM) showed that there 
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Fig I. View from Viking lander I of the surface ol Mars in Chiyse Manilla showing dnits ol fine-grained 
windblown deposits High -resolution images of the drifts show bedding structure, suggesting that the drifts have 
been partly eroded Boulder on the left (‘Big Joe'f is about I x 3 m < Viking Lander I image. M 74301 


were many aggregates of fine-grained particles having the overall 
appearance of miniature popcorn balls (Figure 3). The binding 
force for the particles was strong enough m some aggregates 
that they withstood the vibration of a sonic sifter list'd in 
sieving to determine the particle si/e distribution Further 
inspection, however, showed that some of the aggregates broke 
apart during sieving but then re-formed on the smaller screens, 
suggesting that electrostatic forces were responsible for the 
agglomeration Also, when viewed under a scanning electron 
microscope, the shape of the aggregates would be deformed by 
the electron beam but still remain as aggregates. The possibility 
of electrostatic bonding was reinforced when it was observed 
that the aggregates fell apart when immersed in water where 
electrostatic forces would be neutralized. 

An electrometer was placed on the laboratory apparatus 
during experiments, and an electrical current of 2.4 x 10*^ to 
2.5 x 10'* A was measured in the impact area A comparable 
charge was also observed in wind tunnel simulations (Greeley 
et al . 1^761 for saltatmg sand grains under simulated Martian 
conditions. To determine if the results were artifacts due to the 
experimental apparatus, field measurements of saltating sand 
were made in the Imperial \ alley. ( alifornia. where an electrical 
charge of 1 x 10*^ A was measured for a wind speed slight!) 
above threshold. Samples of sand from the field locality were 
then run in the wind tunnel at I -bar pressure, and a comparable 



fig. 2 Size distributions of particles eroded under simulated 
Martian conditions; velocity of particles was <>1 m . 
atmospheric pressure 5 mbar. Curve on nght shows si/e 
distribution before experiments; alter 20 min of erosion 
the particle size distribution shifted to a finer size (a 
reflection of the kamikaze effect of Sagan et al ||077|i, 
however, as the experiment continued, the size distribution 
shifted to a larger size as the tiny grains began to accrete 
(see f igure 3). 


charge was measured to the field value Similar field results have 
been previously reported by ( till ||04K| I he conclusion is 
thereby drawn that electrostatic charge** occur both on wind- 
blown particles on earth and m simulated Martian conditions 
and that such charges can hold silt and clay -sized particles 
together as sand size aggregates 

Considerable research has been conducted on the effects ol 
electrostatics on small particles; see. for example, summaries 
by /upon ||M0M| and A l> Moore ||973| Although some 
consideration has been given to the effects in the Martian 
environment (Mills. P>77. I den and Vonnegut. 1973), little 
attention has been given to their effects related to geological 
processes Windblown particles in general may be charged in 
many ways, including one or more of the following mecha- 
nisms: (I I contact electrification, resulting from different 

materials coming into contact or from the contact ol the same 
materials, which have different surface properties. (2 )frii tional 
elet trification. resulting from rubbing one surface over another. 
(3 1 piezoclectrii thaws, resulting from pressure. i4)clea\axi 



fig. .V Scanning electron microscope (SI Ml photograph of 
one of larger aggregates from the experiment described in 
figure 2; the ‘popcorn ball* aggregate is composed of 
approximately micron si/e particles held together primarily 
by electrostatic forces (photograph courtesy ol I) k 
Krinsley f 
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fig 4 Current produced hy impact of quartz sand <^250 

pm in diameter ) as a function of atmospheric pressure; 

velocity was held constant 

ein trification, resulting from breakup of certain minerals 
| A I) Moore. 1^73 1 . < 5 ) via thfh ution resulting from frcc/ing 
and thawing of water, and (6) photoelvctrk charging (Pinatti 
and Mascarannas. 1^67) Thus there is ample opportunity for 
windblown particles to become charged. I lectrification may be* 
even more pronounced in the low-density Martian atmosphere 
than under conditions on earth. Paschen's curve shows spark 
breakdown voltage ICobme. 1^41; Llewellyn-Jones. 1 c >5 7 1 as a 
function of atmospheric pressure for COi and earth *ir. 
Although the exact value lor electrical breakdown depends 
upon a number of factors, the curves illustrate that the mini- 
mum values are reached at pressures in the range for the surface 
ot Mars Thus on Mars the atmosphere becomes conducting at 
very low electrical potentials compared to earth 

laboratory experiments were conducted to address the 
question of how electrification of particles could occur under 
Martian conditions and what the effects would be o.i acolian 
processes Three factors were studied: ( I ) the effect of reduced 
atmospheric pressure. (2) the effect of high particle velocities 
associated with strong winds needed to move particles on Mars, 
and (3 1 the effect of partu «c si/e on charge. Figure 4 shows the 
electrical current produced by impact of windblown quart/ 
sand 175-355 /am in diameter on a planar copper plate for a 
range of atmospheric pressures. As pressure is reduced from 
I atm to 4 nibar. the current produced is first quadrupled 
but is very irregular with many ‘spikes’ and then settles at a 
value l 1 : times greater than at I atm. The high initial current is 
thought to result from the fracturing of the grains due to more 
violent collisions which take place between particles when 
cushioning effects of the denser air are reduced at lower pres- 
sure. Initially, the more easily fractured particles were broken, 
generating cleavage electrification. As the test continued, 
frequency of the high current spikes decreased, and only the 
more stable grains remained. During aeolian activity on Mars, 
this would result in higher electrification than on earth, possibly 
with relatively high charges produced in the initial stages of 
wind storms. 


I xpenments were also conducted to determine the etlevt of 
\ locity tm electrification, figure 5 shows that with increasing 
partu le velocity there is a corresponding increase in current, 
again prohahl) resulting front increased energetic collisions 
that resulted in greater ‘cleavage’ land fracturing I electrifica- 
tion. A reversal in electrical sign from positive to negative at the 
higher velocities also occurred 

Changes in electrical sign within the same material were also 
observed as a function of particle si/e I igure h shows current 
as a function of particle si/e t velocity and pressure held 
constant) tor several materials Depending upon the material, 
the sign changed with particle diameter. 

Ihe conclusion drawn from review of the effects ol electro- 
status of particles and from laboratory experiments carried out 
under simulated Martian environments is that < 1 1 electrification 
of windblown particles occurs on earth, (2) under Martian 
conditions, electrification is enhanced because ot reduced 
atmospheric pressure and the high velocities required lor wind 
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I ig 5 C urrent produced by the impact of quartz sand 
('120 n m in diameter) as a function of impact velocity; 
note that there is a change in sign from positive to negative 
at the higher velocities. 
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Fig, 6. Current produced by impact of particles of different 
composition as a function of particle size; note that not only 
is there an effect of composition and size on current but 
that the electrical sign changes as well, 

transport* and (3) the current produced and electrical sign of 
individual particles within a cloud of windblown particles are 
functions of wind velocity (and hence particle velocity)* particle 
composition, and particle size. The aggregates described in the 
first part of this section can be explained as being held together 
primarily by electrostatic Forces, and similar aggregates not only 
may be expected to form on Mars, but their formation is 
probably enhanced. Although such aggregates may be rela- 
tively weak, they are proposed to explain the formation of 
some sandlike structures on Mars, as will be described later. 
First, however, consideration will be given to possible analogs 
on earth, 

Silt-Clay Dunes on Earth 

Dunes and other aeolian deposits composed of silt and/or 
clay have been described in many parts of the world (Bowler 
[1973 ] , Butler [1974] , and others), One of the earliest reports 
of clay dunes is by Coffey [1909], who described dunclike 
ridges many kilometers long, more than 100 m wide, and 10 m 
high on the Texas Gulf Coast, He observed grains of sand 
blowing across the ridges; closer inspection indicated that the 
sand particles were clay pellets that crumbled between the 
fingers but which were strong enough to saltate. The Texas 
clay dunes were studied in more detail by Huffman and Price 
[1949] and, as Coffey had proposed earlier, were found to 
result from a combination of location (proximity to coastal 
brackish lagoons) and climate (dry). Clay crusts were found to 
develop on the dried lagoons where salts had crystallized 


within tlte silt and clay and acted as a cement; as the crusts 
cracked and curled, sharp edges were exposed to (he wind. 
Winds blowing across the dried surface eroded pieces of the 
crust, which, during transport* developed into the sandlike 
‘pellets 1 described by Coffey* The pellets behaved essentially as 
sand as long as they remained dry and were deposited in the 
form of dunes. When wetted, however, the binding forces 
for the pellets were destroyed, and the clay dimes became 
stabilized, 

Similar modes of formation of silt-clay dunes have been 
described in the Colorado Desert [Roth* I960]; Figure 7 
shows a silt-clay dune at (‘lark Dry Lake, California* in which 
the bedding and oilier internal structures are essentially the 
same as ‘normal* sand dunes. Thus even though wetting 
stabilizes silt-clay dunes, the internal structure is preserved. 
Although the process of agglomeration involves primarily 
cementation, and electrostatic bonding of aggregates has not 
been reported on earth (it probably has not been sought), the 
important point is tHl both the external morphology and 
internal bedding structures of dunes and other sandiike forms 
can result from silt and clay aggregates under certain conditions. 

The Mars Model 

It is proposed that some of the dunes, drifts, and other sand- 
like features on Mars are composed of sand size aggregates of 
silt and day particles that are held together by electrostatic 
forces and/or salt cementation, at least during the transporta- 
tion and deposition phases. Several different mechanisms for 
the formation of the aggregates may occur on Mars. 

Derivation From Surface Deposits 

Consideration of the mechanical properties of the fine 
material at both landing sites [II. J. Moore et al„ 1977] 
indicates that it is cohesive and is probably composed of silt 
and day. The term duricrust has been applied to the upper 
layer of these deposits. The binding agent for duricrust could 
be electrostatic* salt cements (salt weathering on Marshas been 
discussed by Matin [ 1974] ), or other interparticle forces, When 
the duricrust is disturbed as, for example* by small mass move- 
ments [Binder et ol», 1977] it displays irregular edges that 
would be susceptible to wind erosion and the formation of silt- 
elay pellets similar to those described by Coffey [1909] for the 
Texas Gulf Coast, Although it may be argued that the pellets 
would be destroyed rapidly under the high winds needed for 
particle threshold on Mars, several opposing factors may 
operate. First, threshold wind speeds needed to move the pellets 
would be lower than for comparable sizes for ‘real* sand 
because, as aggregates* the pellets would be less dense [Greeley 
et al,, 1974a, 1977], In addition, because the pellets would be 
derived from the edges of duricrust or from other aerodynami- 
cally ‘sharp’ places, lower wind velocities could set them into 
motion than could if they were present on fiat beds, Thus lower 
wind velocities would be necessary for pellet movement than 
for real sand, and if the binding strength were high enough, the 
pellets could survive saltation; measurements of the strength, 
however, have not been made. 

A second consideration involves the cascading effect. As 
Bagnold [1941] demonstrated, dynamic threshold speeds 
involving initiation of particle motion through impact by 
material already in motion is considerably lower than static 
threshold speeds. Thus even though the pellets may smash upon 
impact, they could impart enough energy to set other material 
into the airstream which would generate a cascade of material 



Fig. 7. I \ posed section of a silt-clay dune in ( lark Dry Lake. C alifornia, showing bedding structure similar to 
sand dunes. Metal plate is about 7 cm long and indicates scale (photograph by Ron Papton. 1^77) 


moving under relatively low wind speeds. Some of the silt and 
clay derived from the smashed pellets would probably go into 
suspension as a dust cloud. 

Pellet Accretion During Dust Storm Settling 


Aggregates of sand si/e pellets may form from the fine, 
micron si/e particles that occur in dust storms. C ontact and 
frictional electrification of the grains could be expected within 
the dust clouds, given differences in composition, particle si/e. 
and collision.il velocity within the cloud, both positive and 
negative charges would occur. As the dust storm progresses, the 
charges may be cumulative until sufficient charges are reached 


when aggregation begins Once they begin to form, the aggre- 
gates would settle out of the cloud much more quickly than the 
settling rates of the individual dust grams. Formation of the 
aggregates may be an effective mechanism for ending the dust 
storms on Mars As an analog it is well established that dust 
storms on earth are highly electrified | Boning. 1927; kamra. 
|972|. and it has been proposed (Beavers. 1 057 1 that aggre- 
gates of clay minerals were electrostatically attracted and 
deposited as loess in the geological past It is not known, how- 
ever. how large the aggregates might grow 

For the Martian case, as the aggregates settle toward the 
surface, they could be placed into saltation by local winds, or 
gusts, behaving as sand, and be deposited in drifts, dunes, etc. 
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It has been noted [Sagan et aU 19721 that some "variable 
features* such as dark crater streaks and other dark albedo 
patterns begin to develop and are most active following or 
during the waning stages of dust storms. Dark crater streaks are 
generally considered to be eroded zones [Greeley et aL, 1074b; 
Veverka et ah* 19741 where loose particles have been swept 
away to reveal a darker substrate. Dark crater streaks have been 
observed to expand in areal extent to many square kilometers 
on time scales of only weeks, and it is assumed that the wind* 
blown particles eroded are of the size most readily moved, or 
about 160 /4m In diameter. Thus one sequence might be that 
aggregates of sand size grains settle from dust clouds as a thin 
mantle of loose particles: local gusts and turbulence around 
craters and other landforms sweep away the aggregates while 
they are still loosely bonded. With time, in the highly con* 
ductive Martian atmosphere, the charges bonding the aggregate 
would leak* away, and the particles would reassume the 
physical properties of individual particles of silt and clay. 

Accretion During Saltation 

Accretion of pellets may also occur within a saltation cloud. 
Once v iod speeds are sufficiently strong to set either real sand 
or existing aggregates into motion, then the cascading effect 
described above may set other material into motion, including 
fine (silt and clay) particles. Once entrained by the wind* 
electrification could occur within the cloud, which could lead 
to accretion* Even though pellets may be destroyed on Impact 
with the surface, or from interparticle collision* within the 
cloud, aggregates may be constantly reforming in a manner 
similar to that observed during sieving of the laboratory 
samples, The mass of the windblown material, however, could 
hdnw as sand, forming sandlike structures! such as dimes. 

Discussion 

If aggregates of silt and clay form on Mars into sand size 
pellets that are cohesive enough for transport, then several 
interesting effects may occur. First, their formation could 
explain the apparent contradiction between observations 
at the landing sites showing no sand size particles and the 
observations of drifts with internal bedding indicative of sand, 
In addition, some dunes seen from orbit could consist of aggre- 
gates. The mechanism could also explain the relative stability 
of the light streaks in comparison to dark streaks, Light streaks 
are generally considered to be aeolian deposits, whereas most 
dark streaks are thought to be eroded zones [Greeley et ah, 
1974b). Light streaks remain relatively fixed and inactive for 
years, whereas the dark zones shift on time scales of only 
weeks. The white streaks could be deposits of silt and clay that 
were emplaced as aggregates, or pellets that were .sand size 
hence relatively easily moved by the wind, After their emplace- 
ment, however, the binding forces may have been lost (ns the 
case of the clay dunes described by Coffey [19091 and 
Huffman and Price [1949], and the pellets would disaggregate 
in situ, Tire deposit would then reassume its physical properties 
as silt and clay and be very difficult to erode. If the binding 
force were electrostatic, then in the highly conductive Martian 
environment, the charge could be expected to "leak 1 away 
within short time scales. Moreover, the high albedo of the light 
streaks is consistent with fine-grained material such as silt or 
day, 

In addition, aggregates may account for the lower rates of 
erosion on Mars than had been anticipated prior to Viking, 
The high frequency of dust storms observed for many years and 
the high wind speeds needed to move particles led many investi- 


ptors (Sapn [1973] and others to predict very high ndes * 
aeolian erosion on Mars, Viking tmilw* however, indicate rela- 
tively low rates of erosion [Arvidson et nL, 1979]. The low 
rates could be accounted for if many of the saifatmg grains 
arc aggregates, because they would no* iv very effective agents 
of erosion in comparison to whole grains, 

The model proposal here is not considered to apply to all 
Martian aeolian features. Undoubtedly, many of the sand dune* 
observed do involve True* sand pa, ‘ides, not aggregates. The 
mechanism proposed, however, may account for some of the 
observations of the Martian aeolian features. 
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Al burnt- Wind friction threshold speeds (u* t ) for particle move* 
ment (saltation) were determined in a wind tunnel operating at martian 
surface pressure with a 95 percent CO* and 5 percent air atmosphere. The 
relationship between friction speed («*) and l ice-stream velocity Uu) is 
extended to the critical case for Mars of momentum thickness Reynolds 
numbers (Re^) between 425 and 2000. It is determined that lire dynamic 
pressure r uired to initiate saltation is nearly constant for pressures 
between l bar (larth) and 4 mb (Mars) for atmosphere" *if both air and 
CO y, however, the threshold friction speed (u* t ) is alum* J() times higher 
at low pressures than on Barth. f or example, the u* ( (l arth) lor particles 
210 nm m diameter is 022 m s*l and the u n (Mary 5 mb, 200 K) is 22 
mvl 

introduction 

A great deal of attention has been given to lire problem of wind* 
blown particles on Mars, beginning with a paper by Sagan and Pollack 
(1969). Most investigations have involved estimations of wind speeds 
needed to initiate movement of particles of different si/es ami densities 
based on the work originally done by Bagnold (1941) who described the 
physics for windblown sand on Barth, lixirapolatiom were made to Mars 
by substituting the appropriate parameter values for Mars in the equations 
derived by Bagnold For Barth (Sagan and Pollack, 1969, Arvidson, 1972, 
Iverson a al, 1973, and others). Some studies (Chang r/ al, 1968, Adlon 
ct al, 1969) incorporated into the extrapolations data from laboratory 
simulations of some of the conditions on Mars, An extrapolation taking 
interparticle force into account, more appropriate for martian pressures, 
was obtained by I verson ct al ( 1 976)* 

To provide better simulations, an open circuit wind tunnel (Martian 
Surface Wind Tunnel or MARSWIT) was constructed at NASA-Ames 
Research Center (the Facility is described in Greeley ct al, 1977)* lixperl* 
moms from MARSWIT yielded estimates of wind friction speeds for par- 
ticle movement on Mars In which nearly all the important parameters 
were physically simulated or otherwise taken into account. Two of the 
most important parameters for martian simulations are gravity and the 
atmospheric conditions at the surface. For simulations of Mars on Barth 
the primary effect of gravity may be accounted for partly by using pro- 
portionately lower density particles than those expected on Mars. Tire 
atmosphere is more difficult to simulate. Although MARSWIT is capable 
of operating down to about 3 mb pressure, the experiment* presumed 
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by Greeley ct al (1970, 19771 were run not in an atmosphere of martian 
composition (primarily (()*), but rather m ‘larth* air extrapolations 
were made to the martian atmospheric composition by taking into account 
the differences in molecular weight between CO* ami Tartifan However, 
because the dynamics of CO* at low pressures and at the relatively low 
Reynold numbers involved with the martian acolian environment are so 
poorly understood, these extrapolation* were open to question Conse- 
quently, MARSWIT was modified to enable experiments to he run m u 95 
percent CO* atmosphere. The experiments described earlier (Greeley ct al* 
1976) wcrelltiplicaied m the CO* atmosphere to generate a set of threshold 
curved lor particle motion on Mars. The results of these experiments are 
the subject of this report. 

Experimental Procedure and Result* 

The test procedure was us follows; A bed of unifurmwed particles 
was placed m the wind tunnel test section ami the surface was smoothed. 
For each test a 2.5 meter long bed of material was used, as previous tests 
showed (Greeley ct al, 1977) that a bed of this length would give very 
nearly the same results as an mmilelv long hen whereas shorter beds 
resulted m increased saltation threshold values. The tunnel was located 
inside a large vacuum chamber which was evacua ted to about 4 mb pressure 
using ‘Faith' air. The tunnel was then started and the wind speed gradually 
increased until saltation threshold was reached. The tunnel was stopped 
and started several times at each test condition to obtain a number of data 
points. Three separate methods were used to determine sal tat urn threshold; 
(1) A high resolution closed-circuit television system was used to observe 
the particle movement directly; (2) a laser-phot occl system Indicated salta- 
tion by the attenuation of the signal due to particles interrupfing the laser 
beam; and (3) an electrostatic detector measured the current produced by 
saluting particles impinging on the detector element* For most experiments 
ail the systems were used concurrently to detect saltation, with the electro- 
static system being more sensitive for small particles and the television and 
laser systems more sensitive for the larger particles. Many materials were 
tested hut ground walnut shells (specific gravity 1.1 ) were deemed to best 
represent martian particles (Greeley ct al, 1977) and were used to generate 
the data presented here. 

The chamber was then filled with CO* to a pressure of about 80 mb 
giving a 95 percent CO* and $ percent air mixture. Threshold test* were 
run at various pressures as the chamber pressure was again reduced to 4 mb. 
Finally, the chamber was purged witliatr and pumped to 4 mb htr the third 
series ot tests. Threshold speeds were again taken with air as the working 
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Figure 1, Dynamic pressure at the wind tunnel centerline to Initiate parti- 
cle saltation with 5 percent air and 95 percent C0 2 mixture for ground 
walnut shell particles 212 fim in diameter* 


fluid, as the chamber was step-filled to atmospheric pressure to obtain 
threshold speeds for a range of pressures, With this technique it was 
possible to compare a particular bed of material with all conditions nearly 
identical except that the working fluid was changed from air to CO 2 and 
back to air, By comparing the initial air data with the final air data a good 
estimation of the repeatability of the data could be made* Figures 1 and 2 
present typical data sets for a single test bed. Note that the dynamic 
pressure (Ap) is nearly constant for both air and 95 percent CO 2 and with 
chamber pressure (Fig, 1); the velocities (Fig, 2) are different due to the 
differences in properties of air and CC^. This procedure was followed for 
each particle size tested over the range of diameters from 23 /im to 
800 /mi, 

Data Reduction 


The critical measurements in the wind tunnel arc the differential 
pressure (Ap) between the total pressure (p t ) and the static pressure (p s ) at 
the tunnel centerline, and the chamber pressure and temperature, The 
tunnel centerline wind velocities can then be calculated from these measure- 
ments, as well as the friction velocity (u*) for a given bed of materials 
taking into account the boundary layer profile and the surface roughness, 
as discussed in Ivcrscn ctal (1973, 1976). 

Time-averaged boundary layer velocity profiles were determined 
from differential pressure measurements between static pressure and a 
total head pressure tube that was transversed through the boundary layer. 
The tests showed the boundary flow to be essentially steady, two- 
dimensional, constant-property (i.e., maintaining a constant chamber pres- 
sure), and to have a constant free-stream velocity turbulent flow through 
the tunnel, A zero* pressure-gradient flow was achieved by increasing 
the top tunnel wail height with increasing downstream distance (thereby 
expanding the cross sectional area); the zero-pressure-gradient was verified 


by static pressure measurements at several positions in the tunnel. A 
naturally turbulent boundary layer developed for most of the flow condi- 
tions; however, at the lower pressures It was necessary to ‘trip* the boundary 
layer by placing small pebbles on the tunnel floor in the entrance section 
of the tunnel, The test section was more than 25 boundary-layer thickness- 
lengths from the pebbles placed upstream, Insuring a fully developed 
turbulent core region at low pressures. 

The data were reduced taking into account Mach number and slip 
flow effects in the velocity determination. Bach profile was numerically 
curve fitted by means of a multi-plecewise cubic spline technique. 

For the turbulent boundary layer the value of the friction speed u* 
Is dependent on the vaiue of the momentum thickness Reynolds Number 


Re 0 


v 


( 1 ) 


where u w is the free-stream wind-tunnel speed, v is the kinematic viscosity, 
and 0 is the momentum-deficit thickness. The momentum-deficit thickness 
is defined as the thickness of the layer as measured from the wall of an 
external stream (at constant speed u w ) containing a momentum flux equal 
to the loss of momentum flux due to the presence of the wall. U is con- 
venient to define the momentum-deficit thickness 0 in terms of mean 
velocity profiles as 




dy 


( 2 ) 


where 6 is the boundary layer height at which the value of velocity u is 
99 percent that of u*, and y is the vertical height. The momentum thick- 
nesses were determined from integrating the resultant curve fit by a 
numerical quadrature technique. 

An analysis of the zero-pressure-gradient turbulent boundary layer 
for low Reynolds numbers has been made by Coles (1962); from this work 
the surface shear stress r 0 can be determined from the velocity profiles 
obtained in the wind tunnel. This Is accomplished by assuming that the 
velocity obeys the logarithmic law 


in 


M. 
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+ C 
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Figure 2. Wind tunnel centerline velocity to initiate particle saltation 
threshold with air and 95 percent CO^-air mixture for ground walnut 
shells 212 nm in diameter. 
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Figure 3. The ratio of friction velocity (u*) to wind tunnel centerline 
velocity (u*,) us determined by boundary layer surveys and presented as a 
function of Reynolds number (R cx ) based on the distance from the begin* 
ning of the wind tunnel to the test point. 


for the wall region (y < 0.2) excluding the viscous sublayer. In this study, 
k * 0.418 and C * 5.45, as determined by Patel (1965). Converting to 
1°S|0 8* ves 

£ -5.5 log 10 ¥+ 5.45 (4) 

The ratio u^/u^ is well fit by the empirical equation (600 > Re^ > 2000) 

0.702 (log, 0 R cs )‘ hS0 (5) 

as shown by figure 3. In lids equation the relationship is presented as a 
function of £ eX , the Reynolds number based on the distance from the 
beginning of the wind tunnel to the test point. For flows with lower values 
of Rep (Rc 0 < 600) the ratio u^/u^, becomes a constant value of approxi- 
mutely 0.049. This is the maximum value that u yu w cun obtain in a 
turbulent boundary layer at low Reynolds numbers (While, 1979 and Barr, 
1979). For values of Re 0 less than 425 the boundary layer flow is laminar 
and is of no importance in flows of saltating particles. 

Using these relationships the data arc presented as a function of 
particle size versus friction threshold velocity for various pressures (Fig. 4). 
These data are valid for particles larger than about 60 pm, but some adjust- 
ment upward may be necessary for the smaller particles as indicated by the 
following analysis: 

Tests were conducted to determine if the angle of repose, which is a 
function of various interparticle forces, was affected by reducing the 
ambient pressure, Results show that smaller sized particles are affected; for 
example, particles 20 pin in diameter have a much lower angle of repose at 
4 mb than at 1 bar, if the tests were run immediately upon reaching the 
lower pressure. If, however, the particles were allowed to remain at the 
lower pressure for 24 hours before the test was made, the angle of repose 
was the same at 1 bar, Evidently, interparticle air acts as a lubricant as it 
escapes, causing the particles to flow over one another more readily in a 
low-pressure environment. If allowed to remain at low pressure for a 
period of time (many hours), the interparticle air slowly outgassesand no 
longer affects the angle of repose, Ideally, the smaller test particles should 
be allowed to stabilize in the evacuated wind tunnel for a sufficient period 
of time for ail interparticle gas to escape before saltation threshold tests 
commence, This, however, was not possible within the limitations of 
MARSWIT. Thus, the data for the particles less than 60 pm show lower 
thresholds than if the material had been held at low pressure long enough 
for all interparticle gas to escape. (For an analysis of the relationship 
between angle of repose and the threshold velocity see Bagnold, 1941, 
Chapter 7.) Because of this effect the two smaller sizes of particles tested 
showed a decrease in dynamic pressure required for saltation with decreas- 
ing chamber pressure in contrast to the larger sizes tested. This decrease is 
attributed to the lubricating effect of the escaping interparticle gas, and 
correcting for this effect by assuming a constant Ap* the final data curve 
Figure 5 is presented as the best experimental determination of threshold 
friction velocity u* t for particles on Mars, These data are also adjusted for 


representative martian temperatures by the method presented by Iversen 
etaL (1976), 

The particle threshold curves presented here indicate that somewhat 
lower wind speeds are required to move particles on Mars than had been 
predicted by earlier studies. This is due primarily to a more refined method 
of determining the relationship between free stream velocity and friction 
velocity as well as a better simulation using COo as the working fluid. 

Discussion and Comparison with Conditions 
at the Viking Lander Sites 

We have used the results to make a preliminary assessment of the 
ability of winds to move sand sized particles at the Viking Lander Sites. 
To do this, we first relate friction velocity to winds measured by the lander 
meteorology experiment at a height of 1 .6 m above the surface and winds 
at the top of the planetary boundary layer. Both lander sites are charac- 
terized by fairly fiat terrain covered with finegrained material, but littered 
with numerous rocks (Mutch et al, 1976a, b). A reasonable choice for the 
roughness height £ 0 of these surfaces is about 0.1 to 1 cm (Sutton etal, 
1979). Fortuitously, a surface having these values of Z () requires a thres- 
hold wind speed at the top of the boundary layer Vg t that is almost identi- 
cal to the one needed for a surface containing only sand grains (see Fig, 3 
of Pollack et oL 1976a, b). Thus, good estimates of Vg t at the lander 
sites may be made by multiplying the laboratory threshold friction velo- 
cities by the ratio of Vg t /u* appropriate for a sand only case (~50). This 
conversion is indicated by the Case 1 scale of Figure 5. Analysis of the 
Viking meteorology results suggest that the ratio, r, of the wind speeds 
measured by the Viking Meteorology Instrument to those at the top of 
the boundary layer equals about 0,4 (Leovy and Zurek, 1979). 

We are now in a position to compare the velocities measure 
Viking meteorology experiment with threshold values for the op?i n 
grain size, Le„ the grain size easiest to move (—100 fin t), During the 
early portion of the Viking mission (northern hemisphere summer), the 
surface pressures at tiie landers was about 7,5 mb, while close to the 
time of the second global dust storm (southern hemisphere summer) 
the pressure had increased to about 9,5 mb (Ryan et at., 1978). Using the 
results of Figure 5 and a scale factor of 20 to convert friction velocities to 
velocities at the height of the meteorology sensors we estimate that the 
measured wind had to exceed about 30 and 25 m s, respectively, for salta- 
tion to occur. During the summer season in the northern hemisphere, the 
measured wind speed did not exceed 10 ms** at either lander (Hess et ah, 
1976a, b) and hence saltation should not have occurred, in accord with 
lander camera observations (Mutcii et at., 1976a, b). The strongest winds 
measured to date were observed at the first lander site shortly after the 
onset of the second global dust storm, Sustained wind speeds of 18 mf 1 
were recorded, with peak gusts of up to 28 m s'* (Ryan et at., 1978). 
Thus, the winds came quite close to reaching threshold conditions. This 
prediction is in reasonable accord with lander camera results. Movement 
did occur In sand piles placed on top of the lander, where threshold condi- 
tions are lower, but no alteration was observed to occur on the surface, 
except for a few isolated places (Sagan etal, 1977 and Jones et at., 1979). 
Conceivably, during martian years of somewhat enhanced global dust 
storm activity, the saltation threshold is exceeded at the lander sites, 
especially Viking Lander 1, 


P.Jtib 



Figure 4. Friction threshold velocity as a function ui particle size as 
determined in wind tunnel tests using ground walnut shell and 95 percent 
COi-air mixture. 
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figure 5. Particle threshold curves as a function of particle size for 
martian surface pressures and higher, and for representative martian 
temperatures, Case 1 scale is free stream velocity (above boundary 
layer) for winds blowing over a flat surface of crodible grains; Case 2 
surface containing cobbles and small boulders, 
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MARS: PRELIMINARY ESTIMATES OF RATES OF WIND EROSION BASED ON 

LABORATORY SIMULATIONS. R. Greeley and S.H. Williams, Department of Geology and 
Center for Meteorite Studies, Arizona State University, Tempe, AZ 8528 1. 

Ample evidence exists for aeolian erosion and deposition on Mari. Prior to Viking, rela- 
tively high rates of erosion by windblown particles had bee/t predicted (1), based partly on the 
observations of frequent dust storms and the fact that the particles should be highly energetic 
because they would be driven at high velocities by the order-of-magnitude stronger winds requir- 
ed to set them into motion. Furthermore, the decrease in possible atmospheric “cushioning” 
in the low density atmosphere of Mars might enhance the rate of erosion (2). However, one of 
the surprises from the Viking Mission was the apparent low rate of wind erosion evidenced at 
the Viking landing sites (3). A recent estimate (4) suggests a rate on the order of a millimeter 
per million years based oh preservation of surface features and ages derived from impact crater 
statistics. 

In order to determine rates of aeolian erosion on Mats, three general parameters must be 
known: 1) erosive effectiveness of windblown particles under martian conditions, 2) material 
properties of martian rocks and minerals as related to resistance to aeolian erosion, and 3) near- 
surface meteorological conditions on Mars (wind speeds, duration, etc.). The first two general 
parameters can be investigated in wind tunnels and through other laboratory simulations; the 
latter parameter can be obtained from the Viking meteorology experiments, at least for the areas 
around the landing sites.' In this report, we present preliminary data from laboratory simulations 
of aeolian erosion on Mars and combine the results with early Viking meteorology data to 
estimate a rate of erosion. 

The Rotating Arm Mars Erosion Device (RAMED) consists of a holder for various rock 
samples at the end of a rotating arm and a sand hopper from which sand can be dropped into 
the path of the target (fig. 1). Atmospheric pressure and composition, impact particle size and 
quantity, target composition, and impact velocity can all be controlled. The sample is weighed 
before and after each run, the data being plotted as target mass lost per impacting mass, a unit- 
less parameter, as a function of impact velocity and rock type. Quartz sand was used as the 
impacting material. Its size range was 120-180 microns, the size most easily moved by martian 
winds (5). RAMED results for a variety of rock types are shown in fig. 2; as expected, different 
rocks have different resistences to erosion, reflecting grain and crystal size, bonding, composition, 
etc., with rhyolite and obsidian being relatively resistant, and sandstone being relatively suscep- 
tible to erosion. For particles in the 120-180 micron size range impacting at relatively low 
velocities (~10m/s), the mass lost from the target rock is on the order of 10*4 of the impacting 
mass and 10*6 0 f the impacting kinetic energy. An approximate erosion rate can now be calcu- 
lated using the equations presented by Sagan (1). However, several values must be changed to 
reflect Viking data and the laboratory simulations. The frequency of winds of sufficient speed 
to cause particle saltation is several orders of magnitude smaller than that used by Sagan (6). In 
addition, the mass lost from the target per unit impacting KE is on the order of 10*6, rather than 
lO* 4 . A preliminary determination of flux for a saltating cloud was obtained at the MARSWIT 
facility in 1977 (7). When the results are scaled to Mars, the number density of particles becomes 
approximately 10*6 particles cm*3. With these refined values, the rate of erosion is found to be 
about a millimeter per million years and is in general agreement with the estimate of Arvidson 
et al (4), based on a completely different technique. 


WIND EROSION ON MARS 
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The result presented here is preliminary. Additional experiments are in progress to refine 
these data; it is also hoped that additional meteorological data will become available regarding 
wind strengths and frequencies. 
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"STEAM 11 INJECTION OF DUST ON MARS: LABORATORY SIMULATIONS 

Greeley, R. , Center for Meteorite Studies, Department of Geology, Arizona 
State University, Tempe, AZ 8S281 and Leach, R. , University of Santa Clara 
at NAS A- Ames Research Center, Moffett Field, CA 95035 

One of the puzzles regarding the aeolian regime on Mars concerns 
the mechanisms for initiation of the frequent dust storms. Estimated parti- 
cle size In the dust storms is a few microns or less (Pollack et^al. , 1977, 
and others); yet predictions of threshold wind. speeds needed to move particles 
of this size on Mars are in excess of 250 m s~ A (Greeley et_ al. , 1976) for 
winds blowing across a smooth bed of uniform size particles, and it is doubt- 
ful that much, if any, dust is raised in this manner. The optimum particle 
size for movement by winds is about 160pm which has a much lower threshold wind 
friction velocity (*25 to 75 m s" 1 ); once these grains begin saltating they 
could be expected to set smaller grains into motion upon impact. However, 
even these winds seem to be rather infrequent, at least at the Viking landing 
sites. Thus, there may be some mechanism other than wind alone to initiate 
dust movement. Such a mechanism was proposed by Johnson et_ aK (1975) involv- 
ing injection of dust into the atmosphere resulting from desorption in response 
to heating of adsorbed C0 2 on grains. A similar mechanism could occur involving 
water; Huguenin et al^. (1979) have shown that substantial amounts of water 
vapor could be released to the atmosphere and have suggested that some of the 
"blue clouds" observed on Mars are water clouds. 

We have conducted exploratory experiments in the laboratory to study 
the effects of water vaporization on the movement of fine-grained material. 
Experiments were performed in both a bell jar and in a low atmospheric pressure 
wind tunnel using a range of particle compositions and sizes, including: fine 

silt (10pm), talc (13pra), loess (20pm), quartz (23-44pm), walnut shells (61- 
88pm) , and poorly sorted calcite grains (12-1000pm) (walnut shells having been 
used previously in threshold tests, see Greeley et al . , 1977). In bell jar Q 
tests, as the atmospheric pressure was reduced below 10 mb (temperature *v24 C) 
adsorbed water vaporized and ejected particles by one of two processes: 1) vent 

holes and fissureswould develop, followed by a fountainlike spray of particles 
as high as 20 cm above the surface, or 2) violent "eruptions" occurred in a 
boiling fashion (similar to the mechanism described by Johnson et al., 1975). 

The smaller the grain size the more violent the eruption, although some activity 
increased with depth of particle bed, with 20 cm high activity occurring for . 
beds 10 cm deep; some activity was also observed in beds as shallow as 1 mm. 

The amount of adsorbed water also affected activity, with ejection occurring 
with water contents as low as 0.75% by weight. 

Similar effects were observed in tests conducted in the wind tunnel. 
However, it was noted in some experiments that the particle bed remained stable 
until a low velocity wind passed over the surface at which time injection was 
"triggered". Although there is no clear explanation for this effect, the trigr 
gering could be related to slight differences in pressure resulting from the 
wind; 

Figure 1 shows sequential pictures of a silt bed in the wind tunnel; 
as the pressure was reduced, at about 10 mb the surface developed numerous 
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irregular fissures in response to outgassing water vapor. The surface was 
then subjected to wind and the fissured surface began to erode as theyind 
picked up crust-like sections of the silt; the wind speed was *v25 m s , 
substantially lower than "normal" threshold for undisturbed silt. Evidently, 
the fissuring of the silt bed sufficiently roughened the surface to lower the 
threshold speed. On Mars, variations in surface temperature and pressure 
through the critical range for water vaporization may produce "jetting" of 
particles into the atmosphere, or disturb the surface to allow small particles 
to be set into motion by winds otherwise much too slow for "normal" threshold. 

References 


Greeley, R., « al. (1976) Geophys. Res. Letter, ^ < 417-420. 

Greeley, R., et al . (1977) NASA TM-78423, 29p. 

Hugucnin, R. , et al . (1979J NASA Conf. Publ . 2072, 40. (and this issue). 
Johnson, W. , et al. (1975) Icarus 2b_, 441-443. 


Pollack, J., et al . (1977) J. Geophys. Res., 82^, 4479-4496. 



FLOW 

I I I I I I 1 I 

01234567 

cm 

Figure 1. Series of photographs of fine silt particles in a low pressure wind 
tunnel; a) bed prior to experiment; b) fissured surface that occurred as atmos 
pheric pressure in wind tunnel chamber was reduced to about 10 mb; c) surface 
eroded by 25 m s - * wind. 
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Figure 2. Sequential frames showing injection of fine grained particles as a 
result of water vaporiztion at about 10 mb pressure. 
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Material Removal and Production of Fines During Aeolian Erosion of Minerals : 
Mart and Earth 

T. R. McKee and R. Greeley, Geology Department and Center for Meteorite 
Studiea and D. U. Kr ins ley, Geology Department, Arizona State University, 
Tempe, A Z 85281 

Aeolian erosion, or comaunltlon of rock and particle surfaces by 
Windblown materials is an Important process on Earth and certain other planets. 
Laboratory simulations under various conditions are being utilized in a coor- 
dinated program to study a vide range of aeollan processes from geomorphology 
on a megascopic scale to impact erosion on a microscopic scale. Although 
general concepts of aeollan erosion have been developed, the mechanisms In- 
volved in material removal at microscopic scales and subsequent production of 
fine-grained particulate are unknown. The value to planetary geology of an 
analytical Impact erosion model Is its capability to predict trends in aeollan 
erosion when the properties of the materials involved or the existing erosion 
conditions are changed due to new or reinterpreted data. At present there 
have been too few observations of Impact erosion under closely controlled 
conditions to allow development of a complete analytical model. 

In order to develop a general understanding of aeollan erosion, 
single and multiple Impacts of fine-grained material upon rock surfaces are 
being studied at a range of velocities characteristic of Earth and Mar6 (8 - 
107 m/s). Quartz is one of the most common aeollan materials on Earth and 
can be used as a control in a model development program; in addition, its 
material properties and characteristics are well known. In one portion of 
this study, monocrystalline quartz plates and quartz single crystals were 
eroded by quartz sand (125-180 yin dia.) in a sllnger erosion device at known 
impact velocities and impact angles^-. Individual erosion pits were produced 
on rhombohedral faces of quartz single crystals at about a 50° impact angle 
and impact velocities of 16 m/s and 75 m/s. Two types of erosion pits were 
observed: 1) shallow arc shaped pits (Fig. la,b,c) and 2) pits of complex 

morphology (Fig. 2). At 16 m/s shallow pits (Figs. la,b) with an average 
diameter of 15 urn predominate, whereas at 75 m/s the average pit diameter 
is 30 ym and deep, complex pits predominate. The ideal single impact feature 
for a brittle material is characterized by a central deformation zone (impact 
site) surrounded. by a shallow network of radial and circular microcracks^ 

(Fig. 3) which are modified by the impact angle (Fig. 4) such that an incom- 
plete circle is obtained, and by the crystalline structure of the sample such 
that in some cases the circular outline becomes polyhedral (Fig. lc). If the 
erosion pits in Fig. lb and lc were complete, they would both be about 45 ym 
in diameter, producing several flat chips 10-20 ym in diameter and 1-2 ym 
thick. At higher velocities, the central zone of deformation (about 10 Pm) 
would produce a large number of clay-sized particles. The Argonne slinger 
erosion device unfortunately was not designed for collection of the clay-sized 
particles, however this capability will exist with the ASU machine. Alter- 
nately clay-sized particles produced by erosion in a paddle wheel device were 
studied by transmission electron microscopy (TEM) and selected area electron 
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diffraction (SAED). Submicron partidee produced by eroalon at 40 m/a showed 
evidence of plastic deformation, exhibiting asterlsm in SAED patterns (Fig. 5) 
which was not present in the ground quartz used as starting material. 

Plota of quartz erosion rates versus impact angle exhibit maxima 
near 90* which is characteristic of brittle fracture, whereas ductile fracture 
would have a maximum at about 18* (Fig. 6 ) 3 . The surface textures produced 
exhibit angular pits surrounded by irregular, offset ridges (FigT 7, 8a). A 
model for removal of surface material is simply an extension of the complex 
erosion pita produced by interaction of several microcrack systems. Material 
is removed from the surface by the interaction of microcracks radiating from 
a network of adjacent impact sites (circles. Fig. 7) and complex, blocky or 
flat chips 10-30 pm in diameter are produced. At 75 m/s many of the ridge- 
top surfaces appear to be plastically deformed and probably represent remnant 
portions of the central zone of deformation of Impact sites (Fig. 8a, b). 

In summary, it has been shown that aeolian surface textures on 
quartz are produced by interaction of microcracks radiating from individual 
impact sites in a manner characteristic 'of brittle fracture. In addition, a 
mechanism is suggested for the production of fine silt and day-sized material, 
depending upon the Impact velocities involved. At characteristic earth 
velocities silt would probably be the dominant product whereas at martian 
velocities and higher earth velocities both silt and day-sized material would 
be produced. At higher velocities the clay-sized material and some sand and 
rock surfaces would be plastically deformed, altering their bulk response to 
geochemical and biochemical environments. 
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CALIBRATION OF THE MARSWIT TUNNEL FOR DETERMINATION OF PARTICLE THRESHOLD SPEEDS . 

B. R. White, Department of Mechanical Engineering, Unlversltyof California, bavl s , 

CA 95616; R. N. Leach, University of Santa Clara, NASA-Ames Research Center, Moffett 
Field, CA 94035; J. D. Iversen, Department of Aerospace Engineering, Iowa State 
University, Ames, IA 50011; R. Greeley, Department of Geology and Center for Meteorite 
Studies, Arizona State University, Tempe, AZ 85218. 

This abstract describes the procedure and results established In calibration 
of the low pressure open circuit wind tunnel (MARSWIT, described In Greeley et al.., 
1977) for determination of threshold friction speeds. The wind-tunnel test procedure 
was as follows: A bed of uniform-sized particles was placed In the wind-tunnel test 
section and the surface was smoothed. For each test a 2.5 meter long bed of material 
was used, as previous tests showed (Greeley et al . , 1977) that a bed of this length 
would give very nearly the same results as an Infinitely long bed whereas shorter 
beds resulted In Increased saltation threshold values. The tunnel Is located Inside 
a large vacuum chamber which was evacuated to about 4 mb pressure using "Earth" air. 
The tunnel was then started and the wind speed gradually Increased until saltation 
threshold was reached. The tunnel was stopped and started several times at each test 
condition to give a number of data points. Three separate methods were used to 
determine saltation threshold: (1) A high resolution closed-circuit television system 
was used to observe the particle movement directly; (2) A laser-photocell system Indi- 
cated saltation by the attenuation of the signal due to particles Interrupting the 
laser beam; and (3) An electrostatic detector measured the current produced by 
saltating particles impinging on the detector element. For most experiments all the 
systems were used concurrently to detect saltation, with the electrostatic system 
being more sensitive for small particles and the television and laser system more 
sensitive for the larger particles, 

The chamber was then filled with C0« to a pressure of about 80 mb giving a 
95 percent C0 2 and 5 percent air mixture. ^Threshold tests were run at various 
pressures as the chamber pressure was again reduced to 4 mb. Finally, the chamber 
was purged with air and pumped to 4 mb for the third time. Threshold speeds were 
again taken with air as the working fluid, as the chamber was step-filled to atmos- 
pheric pressure. With this technique It was possible to compare a particular bed of 
material with all conditions nearly identical except that the working fluid was 
changed from air to C0« and back to air. By comparing the initial air data with the 
final air data a good estimation of the repeatability of the data could be made. 

The critical measurements in the wind tunnel are the differential pressure (Ap) 
between the total pressure (p t ) and the static pressure (ps ) at the tunnel center- 
line, and the chamber pressure and temperature. The tunnel centerline wind velo- 
cities can then be calculated from these measurements, as well as the friction 
velocity (u*) for a given bed of materials taking into account the boundary layer 
profile and the surface roughness. 

•* 

Time-averaged boundary layer velocity profiles were determined from differen- 
tial pressure measurements between static Pitot tubes and a Pitot tube that was 
transversed through the boundary layer. The tests showed the boundary flow to be 
essentially steady, two-dimensional, constant-property (i.e., maintaining a constant 
chamber pressure), and to have constant free-stream velocity turbulent flow over the 
tunnel wall. A zero-pressure-gradient flow was achieved by increasing the roof 
height with increasing downstream distance (thereby expanding the cross sectional 
area); the zero-pressure-gradient was verified by static pressure measurements at 
several positions in the tunnel. A naturally turbulent boundary layer was developed 
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for most of the flow conditions; however, at the lower pressures It was necessary 
to trip the boundary layer by placing small pebbles on the tunnel floor In the 
entrance section of the tunnel. The test section was more than 25 boundary- layer 
thickness-lengths from the pebbles placed upstream, insuring a fully developed 
turbulent core region at low pressures. 

The data were reduced taking Into account Mach number and slip flow 
effects In the velocity determination. Each profile was numerically curve fitted 
by means of a multi -piecewise cubic spline technique. 

For the turbulent boundary layer the value of the friction speed u* Is 
dependent on the value of the momentum thickness Reynolds number 



v 


when Uoo is the free-stream wind-tunnel speed, v Is the kinematic viscosity, 
and 0 is the momentum-deficit thickness. The momentum-deficit thickness is 
defined as the thickness of the layer as measured from the wall of an external 
stream (at constant speed u» ) containing a momentum flux equal to the loss of 
momentum flux due to the presence of the wall. It is convenient to define the 
momentum-deficit thickness 0 In terms of mean velocity profiles as 


0 ‘ J [c] ( 1 ' £) dV 

where 6 is the boundary layer height at which the value of velocity u is 99 
percent that of u w , and y is the vertical height. The momentum thicknesses 
were determined from integrating the resultant curve fit by a numerical quad- 
rature technique. 

An analysis of the zero-pressure-gradient turbulent boundary layer for 
low Reynolds numbers has been made by Coles (1962); from this work the surface 
shear stress to can be determined from the velocity profiles obtained in the wind 
tunnel. This is accomplished by assuming that the velocity obeys the logarithmic 
1 aw 

“ , 1 In C, 

U*r K V 

for the wall region (y/6<0.2) excluding the viscous sublayer. In this study, 
k s 0.418 and C = 5.45, as determined by Patel (1965). Converting to 1 og 1Q gives 

a . 5.5 log 10 Y * 5 - 45 

The ratio u*/u«, is well fit by the empirical equation 

.1 -1.59 

£ - 0-702 Oog 10 Re x ) 

as shown by Figure 1. In this equation the relationship is presented as a function 
of Re x * the Reynolds number based on the distance from the beginning of the wind 
tunnel to the test point. For flows with lower values of Re§ (Re$ < 600) the 
ratio u*/u«, becomes a constant value of approximately 0.049. This is the maximum 


value that u^/ti* can obtain in a turbulent boundary layer at low Reynolds 
numbers (White, 1979, and Barr, 1979). For values of Re 0 less than 425 
the boundary layer flow Is laminar and Is of no Importance In flows of 
saltatlng particles. 
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WIND EROSION ON MARS: AN ESTIMATE OF THE RATE OF ABRASION. 
S. H. Wflliami and R. Greeley, Department of Geology, Arizona State University, Tempe, Arizona 85281. 

Preliminary estimates of aeolian abrasion rates on Mars have been reported previously (1), 
and the general procedures and approach are presented by Krinsley et al. (2). Earth based 
observations and the global dust storm seen by Mariner 9 led Sagan (3) to predict extremely large 
aeolian abrasion and deflation rates, based on saltation physics and some assumptions of martian 
environmental conditions. This rapid rate is at odds with constraints placed upon it by observa- 
tions of surface features made from Viking images (4,5). Refinement of Sagan’s calculation was 
made based on laboratory studies and Viking data resulting in an aeolian abrasion rate more 
compatible with observations of the surface (1 ). However, values for the parameters in the calcula- 
tion are open to question, and have led to additional work using a different approach, as reported 
here. 

The mass rate of aeolian abrasion is the product of three basic parameters: 

1 ) S^ The susceptibility of a given material to abrasion measured 

in grams target lost per grams impacting mass. 

2) q The mass flux of impacting material on the target, measured 

in grams impacting per square centimeter of target per 
second; q is a function of wind speed and height above a 
surface roughness. 

3) f The frequence of winds of saltation-strength. 

Mass rate of abrasion * * q * f; units: gm cnv^ sec’* 

The first parameter, S A , is measured in the laboratory using abrasion devices that simulate natural 
abrasion under martian conditions (6). The second parameter, q, can be measured for terrestrial 
environmental conditions in the field and in wind tunnels. An extrapolation of q thus determined 
to martian conditions is possible using equations developed by White (7) and can be compared to 
q measured under simulated martian conditions in MARSWIT, a low pressure wind tunnel at 
NASA-Ames Research Center. An estimate of the third parameter can be made from Viking 
meteorological observations. The only other input required for the rate calculation is the relation- 
ship between wind velocity and sand velocity in a saltation cloud. 

The values for all three parameters are currently being refined; however, a first-order 
estimate of the abrasion rate on Mars can be made, based on the following assumptions: 

1) Impactors are quartz sand similar to those used in the abrasion 
devices. 

2) There is sufficient sand on the martian surface to allow a full 
saltation curtain to develop, as in wind tunnel experiments. 

3) Rocks on Mars have similar abrasion susceptibilities as those on 
Earth of assumed similar composition and texture. 

4) The present meteorological conditions as measured the the Viking 
Lander on Mars are similar to those in the past. 
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Values for these parameters can be estimated as follows; 


Parameter 

Typical Value 

Source 

Reference 

S A 

10* 4 

Abrasion machines 

1.8 

q 

0.1 

Earth wind tunnel 
extrapolated to Mars 

7,9 

f A maximum of 10' z > 

more likely 10" 3 

Mass Rate of Abrasion * (10*^) (10* 1 ) (10* 3 ) * 
target density 

Viking observation 

10*® gm cm'2 sec** 

• 10*^ cm sec** 

10 


The results in an abrasion rate of 0.1 centimeters per year (1000 meters per million years). 
Refinement of the measurements of S^, q and f will increase the accuracy of this rate, however, 
most of the values are within order of magnitude and yet results are greatly at variance with the 
observations of the surface that allow an abrasion rate of approximately 1 O' 3 meters per million 
years (4). Thus either the observation is faulty, or the assumptions upon which the calculation 
is based must be in error to explain the large apparent discrepancy. We consider first that the 
impacting material may not be quartz (1 1); although any other “solid" material would also result 
in high abrasion rates. In many areas of Mars, the sediment supply may be limited, so that q is 
much smaller than wind tunnel tests at that wind speed would indicate. Perhaps the assumption 
having the greatest effect is the size of the impacting particle. Sagan et al. (12) suggested that 
sand sized material would quickly self destruct into micron-sized particles due to high energy 
saltation impacts. Observations of atmospheric dust (13,14) and spectral reflectance information 
(15) support the idea that there is a large amount of material in the martian surface environment 
on the order of a micron in diameter, It has been suggested by Greeley (16) that this fine sediment 
forms aggregates bound in part by electrostatic forces. These aggregates may be strong enough to 
saltate and form dunes and drifts, but would disintegrate upon impact rather than abrade surface 
rocks. The efficiency of abrading targets by such aggregates is currently under investigation. 

In summary, sufficient laboratory work has been done on the martian abrasion problem to 
show that there is a major difference between the calculated rate of aeolian abrasion and the 
limits placed on that rate by the observation of surface features at the Viking landing sites. 
Further work is being done to investigate this discrepancy, but it is quite likely that it is due, at 
least in part, to the existance of aggregates of fine particles rather than sand size particles on the 
martian surface. 
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